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The effect of quintessence perturbations on the ISW effect is studied for a mixed dynamical scalar 
field dark energy (DDE) and pressureless perfect fluid dark matter. A new and general methodology 
is developed to track the growth of the perturbations, which uses only the equation of state (EoS) 
parameter wude(z) = Pdde/pdde of the scalar field DDE, and the initial values of the the relative 
entropy perturbation (between the matter and DDE) and the intrinsic entropy perturbation of 
the scalar field DDE as inputs. We also derive a relation between the rest frame sound speed 
Cs,ddb °f an arbitrary DDE component and its EoS u>dde(z). We show that the ISW signal differs 
from that expected in a ACDM cosmology by as much as +20% to -80% for parameterizations of 
iodde consistent with SNIa data, and about ± 20% for parameterizations of iudde consistent with 
SNIa+CMB+BAO data, at 95% confidence. Our results indicate that, at least in principle, the ISW 
effect can be used to phenomenologically distinguish a cosmological constant from DDE. 



I. INTRODUCTION 

It has been known for almost a decade that the uni- 
verse is accelerating [1-4], and considerable evidence has 
accumulated which indicates that the acceleration is due 
to a negative energy component constituting 70% of the 
energy density of the Universe [5-7] . Determining the na- 
ture of this dark energy has become a central challenge 
of cosmology. 

A first step in this direction would be to determine 
whether the dark energy is sourced by a cosmological 
constant (for reviews, see [8]) or a dynamical field (for 
reviews, see [9]). In terms of physical properties, a dy- 
namical field may be distinguished from a cosmologi- 
cal constant by a time varying equation of state (EoS) 
w(z) = p(z)/p(z), and a non-zero (squared) sound speed 
c 2 (z) = 5p{z)/8p{z). As we discuss below, both of these 
quantities could result in an observable signature differ- 
ent from a cosmological constant. 

Both p and p are functions of the field (scalar in the 
case of quint- or k-essence) associated with the dynam- 
ical dark energy, and in principle can be fully deter- 
mined by solving the equations of motion of this field. 
Alternatively, one can parameterize w(z) as a function 
of red-shift without reference to an explicit form of the 
Lagrangian for the field. A representative list of parame- 
terizations is given in [10], and scalar field reconstruction 
from these parameterizations is given in [11, 12]. 

The "sound speed" of a generic fluid is defined as 
c 2 = <5pdde/<5/Odde- The name "sound speed" is a mis- 
nomer for a non-thermal component, and so c 2 (z) can be 
thought of as shorthand for ^pdde/^/Odde- 1 The sound 



1 The thermodynamic squared sound speed is = (8p/Sp)s. For 
an adiabatic system, by definition the entropy S is constant, and 
so the adiabatic sound speed (? a = Sp/Sp. For non-adiabatic 
systems one needs to take into account the non-adiabatic part 
of the pressure perturbation <5p na d, leading to a non-adiabatic 



speed has been a less studied source of insight into the na- 
ture of dynamical dark energy than has the EoS, but its 
effects have been investigated recently in [13-15]. Sound 
speed varies widely among different DDE models. For ex- 
ample, quintessence, with its canonical kinetic term, has 
a constant sound speed equal to unity (in the rest-frame 
of quintessence). On the other hand, k-essence has a non- 
canonical kinetic term, leading to a z-dependent sound 
speed. Other DDE candidates, e.g., dilatons, Chaplygin 
gas, phantoms, and tachyons, each have their own unique 
attributes determining a characteristic sound speed. As 
we show later, the evolution of the sound speed depends 
on the EoS parameter of the dark energy component, as 
well as the dark energy's intrinsic and relative (to the 
matter component) entropy perturbations. 

The dynamical evolution of the field can affect a num- 
ber of physical observables, including the spectrum and 
growth of large scale structure, weak gravitational lens- 
ing, SNIa apparent luminosities, and CMB anisotropics. 
These avenues have been explored in various works [14- 
18]. One particular manifestation occurs in the late-time 
Integrated Sachs- Wolfe (ISW) effect, which measures the 
evolution of the gravitational potential as the Universe 
enters a phase of dark energy domination at z < 2. This 
effect is only significant on large scales (low multipoles), 
since small-scale fluctuations in the gravitational poten- 
tial smooth out along the line of sight. And it is only 
significant at late times since potentials evolve the same 
as the background during matter domination. The late- 
time ISW effect has been detected in cross-correlations 
between CMB temperature anisotropies and surveys of 
large scale structure [19-34]. 

Recent work on the clustering properties of scalar field 
DDE have indicated that, in the context of Einstein's 
general relativity, scalar field dark energy perturbations 
are likely to be anti-correlated to matter perturbations 



component to the sound speed. 
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in the linear regime [35-37]. In [36] it has been shown 
that even if the scalar field is initially homogeneous, it 
eventually acquires a perturbation anti-correlated to the 
matter perturbation as a result of gravitational coupling. 

At low redshifts, very large scale perturbations (on the 
scale of the horizon) are linear and are responsible for 
the late-time ISW signal in the CMB as described above. 
For a ACDM cosmology, the ISW peak in the CMB is 
purely the effect of matter perturbations, of course. The 
question we ask in this paper is how much the ISW signal 
from scalar field DDE linear perturbations differs from 
that due to the ACDM. 

To answer this question we study the role played by 
both the EoS parameter and the (relative and intrinsic) 
entropy perturbations of the dark energy component in 
the evolution of the Newtonian potential. We compare 
effects to the case of a cosmological constant. Our ap- 
proach is different from previous treatments in that we 
apply the equations of linear perturbation theory to the 
case of a generic quintessence component characterized 
by a parameterized EOS w(z;wq,wi) rather than by a 
Lagrangian, allowing for considerable generality in the 
inclusion of a quintessence component of the Universe. 
Our gauge-invariant approach also tracks the evolution 
of entropy perturbations in a consistent manner. 

The outline of this paper is as follows: in the next 
section (II) we describe our mathematical approach 
and demonstrate how it connects with previous treat- 
ments. We summarize two popular parameterizations of 
wdde(z) evolution, and we present a new, "rapid transi- 
tion" parameterization, and then go on to discuss initial 
conditions. Our results can be found in section III, fol- 
lowed by conclusions in IV. 



II. MODELLING THE SYSTEM 

A. Matter and metric perturbations 

We begin with a very general perturbed metric. In 
this Section, we follow the notation found in [9, 38, 39] 
to write 



ds 2 



-(1 + 2A)dt 2 + 2ad i Bdx i dt 
W [(1 + 2V>)<% + 2d l3 E] dx i dx j , 



(1) 



where A,B,ip and E represent metric perturbations and 
a represents the cosmic scale factor. 

We work with a generic cosmic mixture of mat- 
tcr+scalar field DDE, characterized by a total pressure p, 
total density p, averaged velocity potential v, total EoS 
parameter w — p/p and total sound speed c 2 — Sp/Sp. 
We adopt the convention that variables and parameters 
characterizing the total fluid are presented without sub- 
scripts, whereas variables and parameters characteriz- 
ing single components of the fluid are presented with 
component-identifying subscripts. We assume that the 
matter component is a perfect fluid. Several possibili- 
ties exist for the scalar field DDE component - in our 



approach we simply characterize its equation of state pa- 
rameter as wdde (z;wq, wi), where Wo,Wi are arbitrary 
parameters, discussed in Section (HE). 

We now parallel the formalism in [9], presented in the 
longitudinal gauge [40, 41]. This gauge choice corre- 
sponds to a transformation to a frame such that B = E = 
0. In this gauge, the physical gauge- invariant variables 
which characterize the metric perturbations become: 



A 



d_ 
dl 



a 2 (E + B/a) 



A, 



* = -iP + a 2 H{E + B/a) -> -%l> . 
The energy-momentum tensor can be decomposed as 
T$ = -(p + 6p), T° = -(p + p)v, a , 



(2) 
(3) 



T$ = (p + 6p)5$ + I% 



(4) 



where Hp is a tracefree anisotropic stress. Henceforth, 
we assume that the anisotropic stress II is zero. 

The perturbed Einstein equations yield, at linear or- 
der, 



-^<5> + 3H 2 1> + 3H<i> = -AnGSp 



(5) 
(6) 



+ $ = 4irGa(p + p)v (7) 



3$ + 9ff$ 



+ (6H + 6H 2 + -j)* = 4ttG(S p + 36p) (8) 



6p + 3H(6p + 6p) = (p + p)(3$ 



[aHp + p)y]' 

a 4 (p + p) 



A 

H v 

a 

Ua 

a 



Sp 
p + p 



(9) 
(10) 



where a dot, bold or otherwise, denotes a derivative with 
respect to coordinate time t, H = a/ a, and any quan- 
tity preceded by 5 denotes a perturbation in that quan- 
tity. We have also used the relation between <!> and ^ 
in Eq. (5) in the subsequent equations. This relation 
simply reflects our assumption of no anisotropic stress. 
The resulting anisotropic stress-free line element in the 
longitudinal gauge is simply 



ds 2 = -(1 + 2$)dt 2 + a 2 (l- 2$>)dx 2 



(11) 



From this point on we will use the Newtonian potential 
$ to characterize the metric perturbation. 

Since we work with linear perturbation theory, it is 
extremely convenient to transform to Fourier space. The 
Fourier modes evolve independently. For the remainder 
of this paper, we work with individual modes in Fourier 
space. For convenience we suppress the fc-subscripts, and 
ask the reader to keep this in mind. 

We first express the total matter overdensity (S = 
Sp/ p) and the useful combination O = (c 2 — w) 5 en- 
tirely in terms of $ and the background variables using 
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Eqs. (6) and (8), as follows: 
, 2k 2 

o = 



3a 2 H 2 



$ - 2$ 



H 



e 



3H 



2 ■■ k 2 
$ + H(A + 3w)$ + w— $] 



(12) 
(13) 



As is sometimes done (e.g. [14, 35, 41]), we can re- 
express the system in terms of the matter density con- 
trast 5 and the velocity potential v (a scalar field whose 
gradient is the true velocity) as follows: 



k 2 

-me + 3(1 + w)$ - (i + w)— v 



-vH(l 

1 

+- 



3w) 

$ + 



w 



-v 



1 + w 

w . e 

1+w 1+w 



(14) 



(15) 



The growth of the matter overdensity can be shown to 
be governed by the second order equation: 



5(2 - 3w)H 
3(1 



fc2 r k \i 

— wS+ — (1 



w)$ = 



a a" 
w)[& + $(2 - 3w)H] + 3w$ 
i- 2 3H 2 



(16) 



-3HO + 6[ — 7 + — (1 + 9w)] 



a 



Eq. (16) can be deduced from Eq. (10) after substitut- 
ing v from Eq. (7), S from Eq. (12), and using the time 
derivative of Eq. (9). This equation fully characterizes 
evolution in terms of the metric and its perturbation $, 
the overdensity 6, the equation of state parameter w, the 
sound speed c 2 , and the expansion rate H. It can be 
verified that these equations are equivalent to Eq. (30) 
in [41]. 

The pressure perturbation can be decomposed into an 
adiabatic and non-adiabatic component as follows: 



Sp = SpnM + c 2 a Sp 



(17) 



where 5p na d is the non-adiabatic pressure perturbation 
and c 2 = p/p is the adiabatic sound speed of the total 
fluid. 

Finally, following [42], we introduce two gauge- 
invariant entropy perturbation variables. The relative 
entropy perturbation between the matter and the scalar 
field DDE is denoted by the variable S: 



3H(l + lODDE)fifi 

b(z) = - 



l + w 



Sp 



DDE 



PDDE 



Spn 
Pm 



(18) 



The intrinsic entropy perturbation of the scalar field 
DDE is denoted by the variable T: 



r(z) = 



3H(1 + wdde)c 2 dde 



a. DDE 



SpBBE <5pDDE 



PDDE 



PDDE 



(19) 



In terms of these entropy perturbations, the non- 
adiabatic pressure perturbation can be written as: 



SPnad = ^ddeP [(- 



.DDE, 



S 



(20) 



+ (1 



''a.DDE, 



r] 



B. Connection with previous work 

The system of equations (5-10) and the equations (12- 
16) derivative from them are perfectly general, following 
directly from the Einstein equations. We now digress to 
show how what we have done to this point connects with 
previous work. 

Refs. [14] and [35] adopt the approach of starting with 
a set of equations similar to Eq. (14) and Eq. (15), and 
then proceeding with simplifying assumptions regarding 
the behavior of the component (s) of the matter fluid. 

In Eq. (16), if we assume (as is done in [14]) that c 2 = w 
and d(c 2 )/dt = (these assumptions are equivalent to 
setting w, 9 and to zero), we obtain 

5 + 5(2-3w)H + —w5+^r(l + w)$ 

a a z (21) 

= 3(1 + «;)[$ + $(2 -3w)H], 

which coincides with Eq. (375) in [9]. In addition, ignor- 
ing the last two terms of Eq. (12) so that $ is expressed 
by the Poisson equation, and assuming that the evolu- 
tion is largely matter dominated (w ~ 0), one obtains 
the more familiar form 



S + 2HS - 4wGpS = . 



(22) 



This equation is Eq. (66) of [8] and the (unnumbered) 
starting equation in [16]. 

Finally, assuming a late time DDE dominance, defining 
a growth factor g(a) = 5(a) /a in the usual way, and 
rewriting in terms of conformal time and its derivatives 
(denoted by primes), this equation (22) becomes 



5 _ 3w(ffl)OpD E (fl) 
2 2 



+1" 



(23) 



w(a))S! DDE (a) = 0. 



This is just Eq. (1) of [16]. 

In Ref. [15], the authors use Eq. (23) with the met- 
ric perturbation <!> instead of g. Clearly, this equation 
is valid in the longitudinal gauge only if firstly, one 
assumes that $ and 5 are connected through the Poisson 
equation, i.e., $ ~ AnGa^S/k 2 ~ 4:TrGa 3 pg/k 2 and sec- 
ondly, the Universe is largely matter dominated, which 
implies that p ~ a~ 3 . Together, these imply that $ oc <?, 
making Eq. (23) valid for the metric perturbation as well. 

On a side note, Eq. (22), is commonly known as the 
"growth equation" . While this equation is exactly cor- 
rect in the synchronous gauge, it is only approximate in 
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the longitudinal gauge. A detailed study examining the 
accuracy of this equation on different scales in the longi- 
tudinal gauge was performed in [43] , where it was found 
that it can be surprisingly inaccurate on scales larger 
than ~ 0.1ft. Mpc^ 1 . The chief cause for the breakdown 
of this equation on large scales was shown to be the re- 
placement of Eq. (12) by the Poisson equation, and a 
modified growth equation was proposed as a better ap- 
proximation for calculation performed in the longitudinal 
gauge. 



DDE component. Instead, we express the evolution of 
$ in terms of the scalar field DDE EoS wdde ( z 'i wq, Wi), 
and the two independent entropy perturbations S and T. 
Switching to rcdshift z as the time variable, we obtain 
from Eqs. (12), (13) and (20) and Eqs. (20), (33), and 
(34) of [42], 
our central result: 



C. Evolution equation for metric perturbation 



In this paper we avoid making assumptions regard- 
ing the specific nature or behavior of the scalar field 



dz 2 



dS_ 

dz 



d£ 

dz 



1 d$ /3 



1 + z dz V 2 

3 n 



-w — 3c n 



+- 



DDE 



2(1 + *)' 



[-< 



DDE 



2 

S+ (1 



H{zf 



'(l + z) 2 



-a, DDE, 



r] 



l + z 



^3lUDDE - j 



^m C o,DDE 



S + 3 



m m (i 



L a,DDEj 



W 



H 2 
1 



fc 4 (1 + zf ^2(1 + «; DD e) 
ff 4 



l + z 



--(1 + 10)5 + 3 Udde- ., 



9 (l + w) 
1 + w 



k 2 (l + zf 
H 2 



n+,„ ^t? x <? M i k4 ( 1 + z > ( 2 (i + WDDE ) 
(i + '"dde) h--s--t)+ Hi {1 + w) <& 



(24) 



(25) 



(26) 



The different quantities appearing in this system of 
equations are discussed below. It is important to 
note that apart from the evolved variables ((f^S*,! 1 ) 
and 1Z (which is a combination of <f> and d^/dz) all 
the variables appearing in the above equations 
can be expressed as functions of the parameter- 
ized scalar field DDE EoS wdde (z) , as we show be- 
low. For convenience, these variables are listed in Table 
I. 

The quantity 1Z, appearing in the final line of Eq. 26, 
is the gauge-invariant comoving curvature perturbation 
defined as 



n = $ + 



3 (l + w) 



{l + z) Tz 



(27) 



Note that the dynamical Hubble parameter H(z) ap- 
pearing in Eq. (24) is a known function of z: 



H 2 (z) 
H 2 (0) 



= [O m (0)(l + z) s 



(28) 



+ ^DDe(O) £ 



3 J dz 



l + z 



H(0) is the Hubble parameter today, measured to be 72± 
5 km/s/Mpc [44]. The values fi m (0) and £1 D de(0) are 
related by fi m (0) + ^dde(O) — 1 oc the curvature of the 
Universe. Since the first Doppler peak of the CMB offers 
strong evidence that the Universe is flat, we make the 
standard inference that il m (0) + ^dde(O) = 1. 

Moreover, the total EoS parameter w(z) is determined 
by the field EoS, wdde, as follows. In general, w is related 
to the individual EoS parameters Wi of the individual 
components, each with density parameter 0, = Pi/^Pi, 



as 



w(z) 



Qi(z) Wi(z) . 



(29) 



■5 



Hence, for a mixture of perfect fluid matter and scalar 
field DDE, the total EoS parameter w is just 



w(z) = nDDE(z) WDDe(^) 



«dde(0) 



0_ p 3/ o 'dz(i±J£pa) 



H 2 (z) 



wdde(z) • 
(30) 



The adiabatic sound speed of the total fluid <? a {z) is 
computed to be 



2/ \ P 

c a (z) = J= W 



1 + z dw 
3(l+w) ~dz' 



leading to 



1 + w 



(1 + wdde) + 



(1 + z) w[ 



DDE 



WDDE 



(31) 



• (32) 



The adiabatic sound speed of the DDE component is 
given by 



_,„ l w' DDE (z)(l + z) 

"a, DDE ~ »DDE + » — ~ 

6 1 + WDDE 



(33) 



The system of Eqs. (24)-(26) is the central result of 
this paper. The evolution of the metric perturbation (for 
each scale fc) is characterized entirely by a system of lin- 
ear, ordinary differential equations for the evolution of 
Newtonian <&(z), with only the EoS w(wdde(z)), and 
the initial values of the entropy perturbations S and T as 
inputs. These equations present a significant advantage 
of generality in a gauge-invariant formulation. They al- 
low one to track the behavior of different variables such 
as S,v, and, (as explained further below) the dark energy 
sound speed for different classes of dark energy models. 

To the best of our knowledge, the coupled set of equa- 
tions Eqs. (24)-(26) has not been written down in previ- 
ous literature. 

Finally, we derive a connection between the (gauge- 
invariant) dark energy sound speed in the rest frame of 
the dark energy c 2 DDE and the EoS of the dark energy. 
The pressure perturbation of a species "i" in a general 
frame can be related its to the rest-frame speed of sound 
as follows [14]: 

Spi = c 2 Jp t + 3aH(l + Wi)(c 2 s l - c 2 a i )piVi (34) 

Using the definition of T (Eq. (19)), one can deduce 
the following expression for c 2 DDE as follows: 



C S ,DDE ^ c n,DDE 



(35) 



+ 



(1 



a, DDE 



<Sdde + 3aH (1 + wdde) «dde 



Clearly, the evolution of the rest-frame DDE sound- 
speed is linked to the evolution of the EoS parameter 
wbde(z) and the intrinsic entropy perturbation T(z). It 
is straightforward to check that the rest frame sound 
speed of scalar field quintessence is identically unity. 



Note that Eq. (35) is essentially identical to Eq. (8) 
of [45] which was derived in the context of general- 
ized dark matter (subject to a few trivial differences 
in convention- Hu works in conformal time and Thu = 



r us (i 



c a,DDE 



))• 



D. The late-time ISW effect 



The integrated Sachs- Wolfe (ISW) effect is an angu- 
lar variation in the CMB temperature due to a photons 
encountering a time-varying potential well. The relation 
between temperature variation and potential is 



^ AT(n) 



ISW 



2 J d ^ e -frjl(Vo-v)*,v] , 



(36) 

where r{rf) is the optical depth (also called the photon 
opacity), and the integration is along the photon's tra- 
jectory from conformal time r\ r at recombination to the 
present conformal time 770- In this work, it seems safe 
to ignore the small photon opacity t. Then, integration 
of Eq. (36) is trivial, and the resulting ISW relation is 
simply 



= 2($[7to]-$[(7to-7 /r )n,»fc.]). (37) 



ISW 



We have arrived at the ISW relation in Fourier space, 
i.e., as a relation between AT/T and A$ valid for each 
Fourier mode. This relation could be inverse Fourier 
transformed to configuration space, but there is no clear 
theoretical advantage in doing so. 

Thus, to calculate the ISW effect for a given model, 
we are left to evolve the Newtonian potential <& from the 
time of recombination to the present day, using equations 
Eqs. (24)-(26). As inputs, we need the function wdde(z) 
and the initial conditions for the variables $(2) and 
and for the entropy perturbations S(z) and T(z). As we 
have mentioned, the evolution of a potential well does 
not differ from that of a matter-dominated background. 
However, we must integrate from a much earlier time to 
properly include the evolution of the entropy perturba- 
tions S and r. 



E. Parameterizations of the DDE EoS 

To solve Eq. (24), we adopt a phenomenological ap- 
proach and parameterize the input function wude(z), 
rather than derive it from explicit forms of the matter 
action. Several parameterizations have been suggested 
for the scalar field DDE EoS parameter wdde(z)- Many 
of these parameterizations have been shown to faithfully 
mimic the w-behavior of well-known DDE models [46]. 
Such parameterizations are ideally suited for our formal- 
ism, as they allow for a simple deduction of the evolution 
of $ for entire classes of models. In this paper we choose 
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Symbol 


Name 


Defined in 


$ 


metric perturbation 


Eq. 11 


S 


relative entropy perturbation of 
the DDE 


Eq. 18 


r 


intrinsic entropy perturbation 
of the DDE 


Eq. 19 


TZ 


comoving curvature perturba- 
tion 


Eq. 27 


w 


EoS of the total fluid 


Eqs. 29, 30 


^DDE 


parameterized EoS of the DDE 


Eqs. 38, 39, 40 


d 


adiabatic sound speed of the to- 
tal fluid 


Eq. 32 


c a,DDE 


adiabatic sound speed of the 
DDE 


Eq. 33 


r 2 

c s,DDE 


rest frame sound speed of the 
DDE 


Eq. 35 



TABLE I: List of important variables. 



the following two commonly used forms which are suit- 
ably well behaved for large redshifts: 

• The Chevalier-Polarski-Linder (CPL) parameteri- 
zation [47, 48] 

z 

wdde(^) = w + wi —— — - (38) 
(1 + z) 

• The Logarithmic parameterization [49]: 

wdde(z) = w + w 1 \n(l + z) (39) 

In these two parameterizations, wo is today's value w(0). 

Several quintessence models have been proposed where 
the DDE EoS evolves from a steady early value to 
its present value in a rapid transition. The EoS in 
these cases is usually represented by a sigmoidal func- 
tion [46, 50] with four or five parameters. These param- 
eters typically determine the "final" value w(z — 0), the 
"early" value, the point at which the transition occurs, 
and the rapidity with which the transition occurs. As a 
representative of this class of models, we choose a much 
simplified 2-parameter version: 

• The Rapid Transition (RT) parameterization: 

»dde(z) = • (40) 

1 + (W\Z) 

In this parameterization, w is again today's value w(0), 
and \wi\ governs both the location in z of the transi- 
tion and the sharpness of the transition. The evolu- 
tion from initial w = to final wo is half complete at 
zi = l/\wi\. The slope of w(z) at the half-way point is 
w'(zi) = —wq \wi\/2 = —wq/2zi. Thus, a small w\, say 
|wi | < 1, implies an early zi and a relatively weak slope 
(relatively slow transition). On the other hand, a larger 



\wi\ > 1 implies a recent zi and a relatively large slope 
(fast transition). A small w\ (large zi) gives an EoS dif- 
fering little from that of the cosmological constant case 
where w = — 1, whereas a large W\ (small zi) provides a 
considerable difference. 

The functional form for H 2 (z) is given in Eq. (28), and 
that for w(z) in Eq. (30). The same integral exponent 
appears in both equations. Each of the three parameteri- 
zations we use for wdde(z) is conspired to allow a simple 
analytic evaluation of the integral exponent in (28) and 
(30). 

In each of the three parameterizations, wo and wi are 
a priori arbitrary parameters, but later we use the recent 
SNIa and CMB+BAO data to constrain them. For each 
case, we determine numerically the parameter range al- 
lowed at the 95.4% (2a) and 68.3% (la) levels. Then we 
calculate the ISW effect over the resulting allowed ranges. 
In § (III) we will display planes in wq , w\ space with con- 
tours for the allowed 95.4% (2a) and 68.3% (la) regions, 
and contours for the magnitude of the ISW effect. 



F. Initial conditions 

We solve the system of Eqs. (24)-(26) for each of the 
three woDE-pa-rameterizations listed in Eqs. (38)-(40), 
and for a variety of parameter choices (wo ,wi). For each 
parameter choice, we evolve Eq. (24) between last scat- 
tering and now, i.e. a redshift range of z — 1100 to 0, 
eventually filling the {wo, wi}-plane. Then we make an- 
other parameter choice and repeat the procedure, etc. 

Several possibilities exist for choosing the initial (zi ~ 
1100) values of S(z) and T(z). However, in this work 
we focus on the adiabatic choice S(zi) = T(zi) = 0. As 
pointed out in [42] (see also [51] and [52]), it is clear from 
Eqs. (25)-(26) that on supcrhorizon scales (k/aH ^ 1), if 
perturbations are initially adiabatic, then the adiabatic- 
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ity is preserved at all times. Most inflationary mod- 
els ([53], sec e.g. [54] for a review) predict a scale in- 
variant spectrum of adiabatic perturbations, and if the 
DDE perturbations were seeded by inflation, it is reason- 
able to assume that they were initially adiabatic. Ini- 
tially adiabatic modes which are horizon-size today at 
z=0 have been superhorizon (and hence adiabatic) for 
most of the history of the Universe, and hence our choice 
S(zi) = T(zi) = seems well motivated. Of course, it is 
possible to solve Eqs. (24)-(26) to investigate the effect 
of choosing different initial values of the perturbations, 
but we do not do so here. 

For each parameter set, we use the same initial values 
for (<I>fc, $fc), namely, (1,0). As a result of the linearity 
of Eq. (24), the initial choice for is arbitrary; it can- 
not affect the physical results. Put another way, only the 
ratio of the initial and final <i>'s is physical. For the 
choice we have made is reasonable, in that is practi- 
cally frozen during the matter-dominated era when our 
initial conditions arc set. 



III. RESULTS 



Time evolution of variables 



In Figs. (1-2), we show the z-dependence of the vari- 
ables wdde and <f> for the three parameterizations de- 
scribed above. In each case, to visually demonstrate the 
effect, we choose the somewhat extreme parameter val- 
ues wo = — 1 and w\ = 1.4. For the RT parameterization, 
we also examine a more rapid transition having wq = — 1 
and w\ = 3. (In the RT parameterization, wdde is an 
even function of uii, and so the sign of w\ has no mean- 
ing.) All of these choices for {w ,wi} (except for the 
last) are allowed at 95.4% (2a) confidence by the SNIa 
data, but not the CMB+BAO data. 

Fig. (1) shows the behavior of the scalar field DDE EoS 
w(z) for our parameter choices for the different parame- 
terizations. Note that the logarithmic parameterization 
gives physically unrealistic values of w(z) for large z, but 
that is to be expected since this parameterization is un- 
bounded for large z. 

In Fig. (2) find that the late-time dominance of dark 
energy causes the gravitational potential to decrease as 
expected. In a ACDM cosmology, it is known that the 
gravitational potential changes by about 25% between 
last scattering and the present time [15]. This is evident 
from our plot as well. The other curves in Fig. 2 indicate 
that the decrease in $ can be quite different from the 
ACDM scenario. The difference depends on the param- 
eterization scheme chosen, and on the specific values of 
the parameters. 
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FIG. 1: iodde vs z for the three DDE parameterizations with 
(wo, w>i)=(-l,1.4), and also for the RT parameterization with 
(wo, wi) — (-l,+3), and for further comparison, the kja = — 1 
ACDM case. 
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FIG. 2: Potential $ vs z for the three scalar field DDE pa- 
rameterizations with (wq, wi) = (-l,+1.4), and also for the RT 
parameterization with (wo, w i)=(-l,+3), and for ACDM. 



B. Comparison to observations 

We next address the question that motivated this work: 
for a range of "allowed" values of parameters in a given 
parameterization, how different can the ISW effect due 
to scalar field DDE perturbations be from that due to the 
standard ACDM? To answer this question, we compare 
the evolved z = value of <£> from a perturbed scalar 
field DDE scenario to that obtained from evolution in 
ACDM cosmology (where w = — 1 and c 2 = 0). We 
define a "quality" variable Q to characterize the relative 
difference in the ISW effect for scalar field DDE models 
versus ACDM. It is given by 



= 1 



.DDE 

[AT/T] 



(41) 



ISW, ACDM 



For each of the three parameterizations listed in the 
previous section, we evolve the metric perturbations $ 
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and $a over the redshift range z = 1100 to 0. From 
these, we infer the Q values. We fill the {wo, Wi}-plane 
with iso-Q contours, with contour values number-coded 
as shown in Table II. Notice that Q < (Q > 0) means 
that the ISW effect is enhanced (suppressed) for scalar 
field DDE cosmology compared to the ACDM cosmology. 
We will see that typically Q < 0, which implies a larger 
ISW effect in DDE cosmology. 

For each of the three parametcrizations, we then use 
the SNIa standard candle data (ESSENCE+SNLS+HST 
from [55]), to construct a x 2 likelihood indicating which 
values of (wo,Wi) are allowed at the 68.3% (la) and 
95.4% (2a) confidence levels. The x 2 from SNIa is calcu- 
lated as follows: 



N 

J2 [Mobs (Zi) - Mth (zi)Y 
XSN - 



(42) 



where N = 192 is the number of SNIa data points. fj, b s is 
the observed distance modulus, defined as the difference 
between the apparent and absolute magnitude of the su- 
pernova. The a^.i are the errors in the observed distance 
moduli, arising from a variety of sources, and assumed to 
be gaussian and uncorrclated. The theoretical distance 
modulus /ith depends on the model parameters via the 
dimensionless luminosity distance Dl(z): 



H 



fi m (0),ioo,iui) 



(43) 



as follows: 



Mth (z) = 42.38 - 5 log 10 h + 5 log 10 [D L (z)] (44) 

From these results we construct constraint contours 
assuming fl m (0) = 0.28 (best-fit value from WMAP 5- 
year data [5]) and marginalizing over the present day 
Hubble parameter h, following the techniques described 
in [56, 57]. 

Aside from the supernovae standard candle data, one 
can also derive physical constraints from observations 
relating to standard rulers, namely the CMB and the 
Baryon Acoustic Oscillations. Wang and Mukherjee 
showed that together with the baryon density parame- 
ter VL h h 2 , the "CMB shift parameters" [58, 59] defined as 
follows : 

R= yj£l m (0)H o r (z*) , l a = 7rr (z») jr s (z*) (45) 

can be used to set roughly model-independent constraints 
on dark energy models (see also [60]). Here r(z) is the 
comoving distance to redshift z defined as: 



dz 



(46) 



r s (z t ) is the comoving sound horizon at decoupling (red- 
shift z») given by 



poo 



(z) ^3(l + R b /(l + z)) 



dz (47) 



The quantity Rb is the photon-baryon energy-density 
ratio, and its value can be calculated as R b = 
315001M 2 (T CM B/2.7Ky 4 . The redshift at decoupling 
z* is given by the formulas in [61]. 

R can be physically interpreted as a scaled distance to 
recombination, and l a is clearly the angular scale of the 
sound horizon at recombination. 

Following [5], we compute the \ 2 contribution of the 
CMB as 



XCMB ~ VcMB^invVcMB 



(48) 



Here V C mb = P-Pdata, where P is the vector (l a , R, z*) 
and the vector Pdata is formed from the WMAP 5-year 
maximum likelihood values of these quantities [5]. The 
inverse covariance matrix C; nv is provided in [5] . 

A second standard ruler is provided by measurements 
of the Baryon Acoustic Oscillation (BAO) peaks. The 
measured quantity here is the ratio r s (z*) jDy (z), where 
the denominator is the so called "volume distance" de- 
fined in terms of the angular diameter distance Da = 
r (z) I (1 + z) as 



D v (z) 



(l + zfD\(z)z 



1 1/3 



H(z) 



(49) 



So far the BAO peak has been measured at two redshifts, 
z = 0.2 and z = 0.35 [62, 63]. The ratio of the two mea- 
surements of D v (z), i.e., D v (.35) /D v (.2) = 1.812±0.060 
[63], can be used as a model-independent observational 
constraint. In this paper, we calculate the x 2 contribu- 
tion of the BAO measurements as follows: 



Xbao — VB AO C inv V B AO 



(50) 



The vector Vbao = P Pdata, with P = 
(D v (0.32), D v (0.2)) and P data = (0.1980,0.1094), the 
two measured BAO data points [63]. The inverse covari- 
ance matrix is provided in [63] . 

In Figs (3-5) we show the constraint-contours arising 
from both the standard candle (SNIa) and the standard 
ruler (CMB+BAO) data. 



C. Discussion 

The constraint contours allow one to determine the 
range of variation of Q within the regions allowed by the 
observational data. We note that the ISW effects are 
qualitatively similar for the CPL Fig. (3 and the loga- 
rithmic Fig. (4 parameterizations. It is seen that at 2a, 
the Q values range from about +20% to -80% (SNIa 
constraints) and from about 20% to -20% (CMB+BAO) 
constraints. From this one can reasonably conclude that 
our results are valid across a large class of scalar field 
DDE models in which the equation of state gradually 
evolves towards —1 at late times. 
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Number 


Q-range 


1 


Q < -80% 


2 


Q = -80% to 60% 


3 


Q = -60% to - 40% 


4 


Q = -40% to - 20% 


5 


Q = -20% to 0% 


6 


Q = 0% to 20% 



TABLE II: Number-coding for Q values. Negative Q means that the ISW effect is enhanced with scalar field DDE relative to 
ACDM, whereas positive Q means that the ISW effect is relatively suppressed with DDE. 



2 




-1 -0.95 -0.9 -0.85 -0.8 



FIG. 3: Q magnitudes in the {wo, u>i}-plane for the CPL 
parameterization. The number-coding is as in table II. The 
black (thick) lines show the 95.4% (solid) and 68.3% (dotted) 
contours from SNIa data are shown. The blue (thin) lines 
show the corresponding contours from the CMB+BAO data. 




w n 



FIG. 4: Q magnitudes in the {wo, wi}-plane for the loga- 
rithmic parameterization. The number-coding is as in table 
II. The black (thick) lines show the 95.4% (solid) and 68.3% 
(dotted) contours from SNIa data are shown. The blue (thin) 
lines show the corresponding contours from the CMB+BAO 
data. 



The RT parameterization studies the scenario where 
the EoS reaches —1 in a sharp transition, and in this 
case the parameter w\ determines the "sharpness" of the 
jump. Here the data indicate that while sharp (w\ > 2) 
are allowed by the SNIa data, the Q values range between 
20% to -20%. 

The variation of the ISW signal from that expected 
from ACDM for certain ranges of parameters is highly in- 
teresting. In some parts of parameter space for wube(z), 
the ISW signal is enhanced, by as much as 80% for the 
logarithmic parameterization On the other hand, the sig- 
nal is suppressed in some of parameter space. 



IV. CONCLUSIONS 

We have studied the impact of dynamical dark energy 
perturbations on the late-time ISW effect. We devel- 
oped a general gauge-invariant approach for evolving the 
growth of perturbations in the presence of a dynamical 
dark energy. The approach is general in that it uses only 
a parameterization of the scalar field DDE EoS, rather 
than a specific Lagrangian. It also rigorously incorpo- 
rates entropy perturbations in a consistent manner. As 




-0.75 



FIG. 5: Q magnitudes in the {wo, uii}-plane for the RT pa- 
rameterization. The number-coding is as in Table II. The 
RT parameterization is a symmetric function of Wi, and so 
the plot is symmetric about wi = 0; we show only the pos- 
itive half-plane wi > 0. The black (thick) lines show the 
95.4% (solid) and 68.3% (dotted) contours from SNIa data 
are shown. The blue (thin) lines show the corresponding con- 
tours from the CMB+BAO data. 
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an interesting by-product, our formalism allows us to de- 
rive an explicit relationship (Eq. 35) between the scalar 
field DDE rest frame sound speed (? s DDE and ikdde(2), 
which shows that the evolution of these two quantities is 
linked, subject to initial conditions on the entropy per- 
turbations S and T. 

We found that if linear dark energy perturbations with 
adiabatic initial conditions evolve on horizon scales at low 
redshifts, they can enhance or suppress the ISW signal in 
the CMB, depending on the dark energy model used. In 
the case of dark energy models in which the scalar field 
DDE EoS gradually evolves to —1 at late times, the sup- 
pression was as much as 10%, and the enhancement was 
as much as 80% for models allowed by the SNIa data at 
95.4% confidence. Models allowed by the standard ruler 
data (CMB+BAO) were found to enhance the ISW signal 
by 20% or suppress it by 20%. Our treatment can eas- 
ily be extended to other parameterizations of the scalar 
field DDE component and/or different initial conditions 



on the DDE perturbations. 

As an interesting and necessary side issue, we have also 
used the SNIa data to place constraints on a simple 2- 
parameter model in which the DDE EoS can evolve to 
— 1 in a sharp transition. In this instance we have shown 
that while the data do support rapid transitions of the 
EoS parameter, the impact on the ISW effect is restricted 
to within ±20% of the ACDM effect . 
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